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We report the results of dynamical mean field calculations for the metallic Kondo lattice model 
subject to an applied magnetic field. High-quality spectral functions reveal that the picture of rigid, 
hybridized bands, Zeeman-shifted in proportion to the field strength, is qualitatively correct. We 
find evidence of a zero-temperature magnetization plateau, whose onset coincides with the chemical 
potential entering the spin up hybridization gap. The plateau appears at the field scale predicted by 
(static) large- A'' mean field theory and has a magnetization value consistent with that of x — 1 — nc 
spin-polarized heavy holes, where < f is the conduction band filling of the noninteracting system. 
We argue that the emergence of the plateau at low temperature marks the onset of quasiparticle 
coherence. 



I. INTRODUCTION 

The paramagne tic, metallic ground state of heavy 
fermion compounders! can be interpreted as a Fermi liq- 
uid with low coherence temperature and large effective 
mass.l^ The origin of the coherence temperature lies in 
the nature of the quasiparticles, which emerge as a co- 
herent superposition of the Kondo screening clouds of 
the individual magnetic impurities. Being a Fermi liq- 
uid, the heavy fermion ground state is well captured by 
large- iV mean-field theories that predict a renormalized 
band structure; at this level of approximation, the key 
properties — such as the effective mass and the Fermi- 
surface topology — are reproduced. 

Mean-field approaches that are static in time, however, 
do not properly incorporate Kondo screening and hence 
do not account for the many-body nature of the quasipar- 
ticles. Such approaches become questionable when the 
system is subject to perturbations, such as large temper- 
atures or strong magnetic fields, that have the potential 
to destroy Kondo screening — hence, the motivation to 
consider dynamical mean field theories (DMFT), where 
the Kondo effect is built in. 

The problem of a Kondo lattice insulator (symmet- 
ric conduction band at half-filling) in a magnetic field 
has previously been treated approximatelyusing 

DMFTS! 

and exactly using quantum Monte Carlo.'^ In that spe- 
cial case, the only important low-energy scale is the indi- 
rect gap separating upper and lower quasiparticle bands. 
In the metallic case (e.g., filling less than half), there is 
an additional, smaller energy scale given by the separa- 
tion between the chemical potential and the top of the 
lower band.—' This scale is generically small regardless of 
the conduction band filling and controls the Fermi liquid 
properties. As we shall see, it also sets the magnetic field 
scale for the onset of metamagnetic features. 

The static mean field scenario for application of a 
magnetic field at zero temperature is as follows. As 
the field strength is increased, the spin up quasiparti- 
cle bands descend. The spin up Fermi surface shrinks 
to a point and disappears as the lower band drops 
below the chemical potential.'^EIl2] The resulting half- 



metallic stat^iilhas mean field parameters that are locked 
at the values obtained before the disappearance of the 
spin up Fermi surfa ce an d are completely insensitive to 
changes in the field .1^^ As a result, the physical mag- 
netization —dJ-/dB (where T is the free energy) has 
a constant slope proportional to the difference between 
the Lande g-factors for the two species. If the c- and 
/-electrons couple identically to the applied field, the 
system is predicted to exhibit a magnetization plateau 
at a value that depends only on the conduction band 
filling ric. Note that the dimensionless magnetization 
M = (uc,! — nc,x)-l- (n/j — n/,|) always shows the plateau 
behaviour, irrespective of the values of gc and gf. This 
is related to the fact that, in the locked state, the sys- 
tem behaves as a gas of a; = 1 — fully spin-polarized, 
heavy-quasiparticle holes. 

Our DMFT results verify this scenario. We find that 
the magnetization profile begins to develop an infiection 
at temperatures on the order of the predicted coherence 
temperature. This feature resembles a plateau smeared 
by temperature, and we have verified that it sharpens 
as the temperature is lowered. Moreover, the appar- 
ent height of the plateau is consistent with the predicted 
magnetization value. An analysis of the evolution of the 
spectral functions with applied field shows that the end- 
points of the plateau do coincide with the chemical po- 
tential entering and leaving the hybridization gap. 

II. MODEL AND METHODS 

We consider the Kondo lattice model (KLM) on a 
square lattice in an external magnetic field B: 

' (1) 
Ho = Y^ieu - M)4,sCk.. - gi^BBY^iSl + Slj. 

k,s i 

Here, ^ creates a conduction electron in an extended 
orbital with wavevector k and spin projection s = Tj i 
along an axis of quantization chosen parallel to the 
applied field. The tight-binding dispersion relation is 
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Ek = —2t{coskx + cos ky). At each lattice site i, a lo- 
cal spin-1/2 degree of freedom S/ is coupled via J to the 
c-electron spin density S? = 5 X^s s' '^i s'''s,s'Ci,s' (repre- 
sented with the aid of the Pauli spin matrices a). An 
analogous expression can be written for S/ using the lo- 
calized orbital creation operators Z/^. Since the KLM 
forbids charge fluctuations on the /-orbitals, this repre- 
sentation demands that we impose a strict constraint of 
one electron per /-orbital. 

The DMFT approximation neglects spatial fluctua- 
tions and thereby omits the k-dependence of the self- 
energy: i.e., Ss(k, ia;,„) Hs{iujm).^'^ Being local, the 
self-energy is that of a single impurity in a bath described 
by the free Hamilton operator Tio = s ^k.scjk s'^k.s — 
gf^BBSf^. The corresponding local Green bath function 
is given by 



Gri^^m) 

GiHliOm) 



(2) 



using g'f{iuj,n) = -/o^dre''^"^(co,,.(T)4 ,)Ho and the 
equivalent definition for Q^^ . The prerequisite for the 
implementation of DMFT is the ability to solve, for a 
given bath, the Kondo model. 



TL = Ho + JSq • Sq 



/ 



(3) 



In our calculations, we have opted for the Hirsch-Fye 
approach. Following Ref.[T4l we expand the Hilbert space 
to permit charge fluctuations on the /-sites and replace 
the Kondo term by 



h.c. 



+ ^(/aT4T-l/2)(/^,i/o.x-l/2). (4) 

At the expense of two discrete Hubbard-Stratonovich 
fields, we can decompose the (perfect square) hybridiza- 
tion and Hubbard terms, and readily implement the 
above interaction within the framework of the Hirsch- 
Fye algorithm.lSl In the limit U ^ 00, charge fluctu- 
ations on the /-orbitals are suppressed and the squared 
hybridization reduces to the desired Kondo term. The ef- 
flcicncy of this approach lies in the fact that C/(/q |/o,t ^ 

1/2) (/q |/o.i — 1/2) is a conserved quantity^"* for the con- 
sidered class of Ti.0 , which does not include hybridization 
terms between the bath and impurity orbital. Typically, 
f3U ~ 30 suflices to ensure that double occupancy drops 
down to (/d.|/o,T/d.i/oa) = 0-0005 ± 0.0005 and that 

{J2s fo.sfo,s) — 1 to the same precision. Since charge 
fluctuations on the /-sites are suppressed by the Hub- 
bard interaction, the local Green function, as obtained 
from the Hirsch-Fye alogorithm, is diagonal in the large- 
U limit and is given by 



GsiiuJm) 
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(5) 



Self-consistency requires that the local Green func- 
tion, as determined from the effective impurity problem, 
matches that of the lattice: 



1 



1 



(6) 



Here, Go.s(k, iujm) is the free lattice Green function [cor- 
responding to Hq in Eq. 

Hence, for a given bath Green function Qsii^m), we 
extract the self-energy from the impurity solver and, ex- 
ploiting Eq. ^ , recompute the bath Green function with 



^VGoi(k, 



(7) 



The procedure is repeated until convergence is achieved. 
As noted above, in the limit [/ — 00 the local Green 
function, self-energy, and bath Green function arc all di- 
agonal. Accordingly, the self-consistency [see Eq. (|6|] is 
equivalent to two scalar equations: 



k 0,s 



"^(k,iw™) - (iwm) 



1, 



(8) 
(9) 



The last equality follows from the fact that GQ{(k, iujm) 
has no k dependence. 

Having determined the self-energy, one can readily 
compute the single particle Green function Gf{}i^ iujm) 
and extract the spectral function A(k, uj ) with a stochas- 
tic analytical continuation method.'SEH 



III. MEAN FIELD PICTURE 

Free electrons moving via nearest-neighbour hopping 
on the square lattice have a density of states (DOS) 

poiuj) = ^^(Vi - {io/uy)eim' - lu^), (lo) 

written in closed form in terms of an elliptic integral of 
the flrst kind, JC{x). The function has support in the 
region —it < u < At and thus a bandwidth W = 8t. 
It is smooth and continous, except at the step-like band 
edges and at w = (a van Hove singularity), where it 
diverges as /9o(u;) ^ SttI* T^f ■ [Numerical evaluation 



of Eq. ( 10 1 can be carried out very efficiently using series 
expansion; see Appendix [A|] 

At the mean-field level, the interacting Kondo lat- 
tice system is modelled by a bilinear effective Hamil- 
tonian that includes both c- and /-electrons. This 
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FIG. 1: The upper and lower hybridized bands are plot- 
ted over a trajectory bounding one octant of the Brillouin 
zone. The grey line-width is proportional to the c-electron 
spectral weight. The indirect gap 2 A and the optical gap 2V 
are indicated. The corresponding c-electron density of states 
pc{i-o) is plotted to the right on the same energy scale. In 
magnification, the region around the gap is illustrated with 
the filled Fermi sea shaded in grey, revealing a small amount 
of headroom et <Si 2A in the lower band. In the lower two 
panels, the small, noninteracting Fermi surface is contrasted 
with the large interacting one. 



can be obtained from the saddle-point of the Hubbard- 
Stratonovich decomposition of Eq. (|4| : 



-V 



-V \ 



(11) 

In the usual way, hybridization between the c-electrons 
and the dispersionless band of /-electrons leads to a 
renormalized, quasiparticle dispersion 




rir. 



FIG. 2: (Upper panel) As a function of the bare filling ric, 
the inverse quasiparticle DOS l/p(/x/) and the indirect gap 
2A are compared to the energy scales defined by the Kondo 
temperature Tk (with fcs = 1) and the width et of unfilled 
states at the top of the lower band. (Lower panel) The dashed 
lines at 2et and 2(et -1- A) mark the plateau endpoints that 
arise from a purely rigid shifting of the bands (in units where 
fiBQ = 1). The shaded areas between B' and B" indicate the 
extent of the locked region as calculated using the fully self- 
consistent mean field equations. The discrepancy is a result of 
higher order changes in the mean field parameters. We expect 
that the inclusion of V phase fiuctuations restricts plateau 
formation to between B' and B' + A' (roughly) and leads 
to a slow disintegration of the heavy fermion state between 
B' + A' and B" . The vertical dashed (red) lines between B' 
and B' + A' at ric — 0.65 and Uc = 0.85 appear again as 
horizontal dashed (red) lines in Fig. [g] 
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ek - sgfiBB ± V(ek - 6)2 + 41^2 



(12) 



and a quasiparticle DOS p^uj) — pc{^) + Pfi'^)i where 



Pc(w) = po{uj - V^/uj + h) 



(13) 



and Pf{i^) — {V'^ / uj"^) pc{io) . Here b — p,c — Pf is the 
chemical energy to transmute the c- and /-electron char- 
acter of a particle, and is a hybridization energy de- 
termined self-consistently via V ~ J^si^k sf^,^) ■ The 
spectral weight vanishes outside the lower [loi , ^2] and 
upper [ti;3,a;4] hybridized bands, where LUi {i = 1,2,3,4) 
denote the four ordered roots of ±4i — uj — V"^ j uj + 



The mean field equations can be written compactly as 



E 



I 



-J/2uj 



(14) 



where the integral is taken over the disjoint interval / = 
[wi, aj2] U [0^3, 0^4], and pf^s is shorthand for pf ±gpBB /2. 
The first two equations fix the c- and /-electron occupa- 
tion (to Tie and 1, respectively), and the third enforces 
the self-consistency condition on V. We have written 
J' = (const.) X J to allow for the fact that different mean 
field decompositions lead to a different numerical prefac- 
tor. In each spin channel, 2A = tJa — 0^2 is the small- 
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FIG. 3; The magnetization of conduction rric and local mo- 
ments ruf, computed via DMFT, are plotted as a function of 
magnetic field at various temperatures. In the upper panel, 
the lines running through the /-electron data points are fit 
to the forms tanh(S"=f^p^isB/3). At f3t = 5, 23"" = 1.03 and 
at I3t = 10, 25"^^ — 0.65. In the lower panel, the magnetiza- 
tion profile is marked with three regions. A sketch of the spin 
resolved density of states as deduced from the single particle 
spectral function (See Fig. [5| is given in each region. 



est indir ect ga p, and 2V is the threshold for optical ex- 
citations .'^^'^ Figure [l] illustrates the zero-temperature, 
zero-magnetic-ficld solution appropriate for band filling 
< ric < 1. An additional energy scale et = UJ2 — fJ-f, 
representing the headroom at the top of the lower band, 
is indicated. 

An artifact of the mean field treatment is that the hy- 
bridization matrix element has an anomolous expecta- 
tion value that vanishes with heating at a second-order 
phase transition. (The true Kondo physics is that of a 
crossover.) There is a critical temperature Tc such that 



as T ^ Tc from below, V^/i/ — > and b fic- As 
the quasiparticles disintegrate into their separate c- and 
/-character consituents, the corresponding densities of 
states return to their free values: Pc{^) Pq{'-^ + A*c) 
and Pf{uj) S{iL!). In this hmit, Eq. ([Til) reduces to 



2 
J 



duj 



Po(w) 



tanh[(ti; — pc) 721^ 



LU - He 



UJ - Pc 



Po(Mc)log^f— ^, (15) 
4t -f- Pe 



where Pc{nc) takes its noninteracting value, determined 
implicitly by 



2 / po{Lj)f{u; - pc) = ric. 
I -At 



(16) 



The critical temperature scales as ~ aW , where 
the constant of proportionality is a function of the band 
filling alone. (For a flat DOS, Tc = 0.567^^ x aW .) 
The small parameter a = e^^/^''°^'^''^ renormalizes the 
bandwidth down to the Kondo scale; hence it is natural 
to identify a Kondo temperature Tk = T^. On the square 
lattice, the value of Tk is always nonzero for Uc > Q. Tk 
is compared to the important zero-temperature energy 
scales in the upper panel of Fig. |2] A value J /t = 1. 631 
was chosen so that dM / dB\^_^ = p{Pf) matches the 
DMFT results for J/t = 1.6 and = 0.85. (We set 
gpB — 1 for the remainder of this section.) 

At zero temperature, where the hybridization is at its 
strongest, the Fermi function in Eq. ( [l4| ) cuts off the 
integration at /i/,s. There are three possibilities: (i) cui < 
Pf^s < W2; (11) wi < pf^i < UJ2 and UJ2 < Pf,] < W3; (HI) 
LOi < pf^^ < UJ2 and Wa < pf^^ < LU4. In cases I and 
III, both spin up and spin down quasiparticles have a 
Fermi surface. In case II, the chemical potential lies in 
the spin up hybridization gap. As a consequence, the 
upper limit of integration in the spin up channel is lu2, 
and no longer enters the equations. The magnetic 
field enters only indirectly through /i/,x, and thus the 
system becomes insensitive to changes in magnetic field. 
In particular, this leads to a plateau in the magnetization. 
Note that pf,i = const, and hence pf,^ = const, -t- B, 
which means that the upper band descends at twice the 
rate it does in case I, where pf^i^ ^ Pf\B=o + B/2. The 
plateau terminates when pf^^ reaches the lower edge of 
the upper hydridized band. 

The leading edge of the plateau can be determined 
by solving for the field B' at which | = 0^2 and 
Pf,l = UJ2 — B' . Since the mean field parameters are 
locked beyond B' , the far edge is simply B" = B' + 2A', 
where 2A' is the size of the indirect hybridization gap 
at B = B' (and its size throughout B' < B < B"). 
This argument for the location of the far edge is flawed 
in one respect: when B > B" , the system is expected 
to enter region III, but the mean field equations in that 
region turn out to be pathological {V begins to increase 
rapidly with B). Indeed, energetic considerations tell us 
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that for B = B"[l - 0{a'^)] the V = normal state 
becomes energetically favourable.^ Nonetheless, this ap- 
parent first-order collapse to the normal state cannot 
survive the inclusion of phase fluctuations in V, which 
smooths out the destruction of the heavy fermion state 
over a width A. Hence, we expect to find a plateau in 
B' < B < B' + A' and a crossover to the normal state in 
B' + A' < B < B" = B' + 2A'. The predicted extent of 
the plateau is shown in the lower panel of Fig. |2] 

The mean field picture describes x = 1 — Uc heavy 
holes, which spin-align with the applied field and become 
completely polarized at B' . Along the plateau, the value 
of the magnetization is M = (ricT + ^/t) ~ ("-cT + ^/t) ~ 
x[\ — 0(a^)]. This follows because ric] — Uf^ ~ 1 and 
rici — '"^fi — 1 ~ 2;/2 when the chemical potential sits 
in the spin up hybridization gap. In the Kondo limit 
(i/ ^ W), a is exponentially small and corrections of 
order are completely negligible. Hence, the plateau 
has a height of almost exactly M — x. 



IV. DMFT RESULTS 

The mean-field picture of the metamagnetic transition 
and the associated change in the Fermi surface topology 
is a direct consequence of coherence. To confirm this, we 
have carried out temperature and magnetic field DMFT 
scans at J/t =1.6 and electron density ric = 0.85 (see 
Fig. [3]). This choice of parameters sets the single im- 
purity Kondo scale to Tfc/t ~ 0.09.123 At temperature 
scales larger than the Kondo scale, one expects the local 
moments to be essentially free with magnetization 
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A(k,a;), qubB^O 



(17) 



where 5®*^ allows for a renormalization of the impurity 
moment. At fit = 5, the DMFT data (see top panel of 
Fig.pl for nif follows this form with S'^^ ~ 1/2 to consid- 
erable accuracy, thereby showing that the system is in the 
high temperature local moment regime. The polariza- 
tion of the conduction electrons is opposed to the applied 
magnetic field because the antiferromagnetic Kondo in- 
teraction generates a negative effective magnetic field .'^'iSI 
In the vicinity of the Kondo temperature, at (3t = 10, 
the magnetization curve to/ can roughly be accounted 
for with the free moment form of Eq. (17 1, albeit with 
geS ^ This reduction is consistent with the onset 

of Kondo screening. At temperatures (3t > 30, which we 
argue are well below the coherence temperature, three 
distinct regions denoted by I, II, and III are apparent 
in the data, as shown in Fig. [3] Our understanding of 
those features relies on the single particle spectral func- 
tion below the coherence temperature and the associated 
change in the Fermi surface topology as a function of the 
magnetic field. 

In zero magnetic field and at our lowest temperature, 
Pt — 40 [see Fig. 4](c)], the single particle spectral func- 
tion follows the hybridized band picture of the mean- 
field calculation (cf. Fig.[T|), refiecting the coherent heavy 



uj/t 



uj/t 



Lo/t 



0.001 



0.01 



0.1 



6- 
4- 
2- 
0- 
-2- 
-4- 
-6- 

6- 
4- 
2- 
0- 
-2- 
-4- 
-6- 



1 10 

(a) /3i = 5 



(b) [it = 10 



4 



(c) (it - 




(7r,7r) 



FIG. 4: The c-electron spectral function v4(k, uj) in zero field 
is plotted for three progressively lower temperatures. The 
legend at the top indicates the false color values. 



fermion metallic state with Luttinger volume including 
both conduction and localized electrons and large effec- 
tive mass. Introducing a magnetic field leads to com- 
peting effects. One possibile outcome is that Kondo 
screening is completely suppressed — thus triggering the 
breakdown of the heavy-quasiparticle, the heavy fermion 
coherent state, and the associated large Fermi surface. 
Another is that partial Kondo screening persists and 
that the heavy fermion metallic state survives in some 
form. Fig. [s] s upports the latter. In particular, in re- 
gion I [Fig. [5|^a),(b)] the two-fold spin degeneracy of 
the spectral function is lifted and the spin down (up) 
band shifts to higher (lower) energies, thereby producing 
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(a) Ai{k,uj), g/isB^ 0.0125 



(b) ^t(M), 3^^35 = 0.0125 
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(c) A|(k,c^),5A*BS = 0.05 




(e) A|(M),5a*bS = 0.12 
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FIG. 5: Low temperature, /3t = 40, spectral function in the up (left) and down (right) spin sector as a function of magnetic 
field. We use the same scale as in Fig. |4] 



hole-like Fermi surfaces of different size centered around 
k — (tt, tt). [The fact that the spin down (up) band shifts 
up (down) in energy and yet rric < is a consequence of 
some subtle rearrangement of spectral weight.] The on- 
set of region II [Fig.[5jc),(d)] is marked by the vanishing 
of the spin up Fermi surface: the chemical potential lies 



within the hybridization gap, and the lower spin up band 
is completely filled. This sudden change in the Fermi sur- 
face topology is the origin of the cusp-like feature sepa- 
rating regions I and II. Note that the magnetization data 
of Fig. |3] supports the sharpening of the cusp as the tem- 
perature is lowered. Region III [Fig. |5je),(f)] is again 
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0.15 



FIG. 6: Total magnetization, M — nic + nif, as function 
of temperature and band filling. The plateaux at a:: = 1 — 
ric, predicted by the mean field theory, are indicated with 
horizontal dashed (red) lines and correspond to the vertical 
dashed (red) lines in the lower panel of Fig. [2] 



characterized by a change in the topology of the Fermi 
surface. Here the spin up valence band drops below the 
chemical potential, thereby forming an electron-like spin 
down Fermi surface centred around the k = (0,0) point. 
Across regions I, II, and III, the spin down Fermi sur- 
face evolves continuously with a growing hole-like Fermi 
surface centered around k = (tt, tt). A sketch of the spin- 
resolved density of states in the three regions is given in 
Fig. [3] At very high magnetic fields [Fig. |5jg),(h)] the 
spectral function progressively tends towards that of a 
Zeeman split cosine band. 

The onset of the plateau-like feature in the magneti- 
zation curve as a function of temperature can be used as 
measure of the coherence temperature. At J/t = 1.6 and 
Tic = 0.85, Fig. |3] provides a rough estimate of this scale: 
Tcoh/t — 1/30. Given Tcoh, it is interesting to analyze the 
temperature dependence of the single particle spectral 
function (see Fig. [3]). Already at temperatures T/t = 1/5 
and T/t = 1/10 above T^/i ^ 0.09, features of the hy- 
bridized bands are apparent. At those temperatures the 
magnetization curve is featureless since the coherence or 
hybridization gap is not formed. At T/t = 1/40 < T^oh 
(Fig.|4|;) the coherence gap is well formed and the plateau 
feature in the magnetization curve is apparent. 

As argued in Sec |III| the height of the magnetization 
plateau, M = 1 — ric, follows from the Luttinger sum 
rule. Fig. [6] plots the total magnetization M — rric + ruf 
at two different conduction band fillings, Uc = 0.85 and 
ric = 0.65. As the temperature decreases, a plateau-like 
feature precisely at 1 — ric emerges. For comparison, we 
have plotted the mean-field prediction as obtained from 
Fig. |2] Note that the initial slope in the magnetization 
curve is related to the effective mass, which is inversely 
proportional to the coherence temperature. Hence, at 



low band fillings, lower temperature values are required 
to reveal the plateau feature. 



V. CONCLUSIONS 

The interaction of conducting electrons with local im- 
purities often gives rise to complicated nonlinearities in 
magnetic response (often described as metamagnetism) . 
Below the coherence temperature, we find that the mag- 
netization profile of the Kondo lattice system shows a 
plateau, whose onset is linked to the vanishing of the 
Fermi surface in one spin sector. The relevant energy 
scale corresponds to the headroom above the chemical 
potential in the lower hybridized band. As a consequence 
of the Luttinger sum rule, the total magnetization is 
locked to the value M = 1 — ric throughout the plateau. 
The plateau survives up to magnetic fields at which the 
quasiparticles begin to break apart. At high fields, the 
system is once again characterized by Fermi surfaces in 
both spin sectors, which smoothly evolve with increasing 
field back to those expected for the Zeeman-split bare 
conduction band. 

We have reached these conclusions on the basis of dy- 
namical and static mean field calculations, both of which 
are consistent with the scenario we have outlined. There 
is one important disagreement between the two theories. 
At large fields, the DMFT shows that once the chemi- 
cal has traversed the gap, it enters the upper hybridized 
band, and the system again becomes a heavy metal with 
two Fermi surfaces. The behavior of the static mean 
field equations after the reentry into the upper band is 
pathological. Energy considerations suggest a first-order 
collapse back to the normal state. There is no evidence of 
this within DMFT, and we expect that phase fluctuations 
of the hybridization field are responsible for washing out 
this false transition. 

In closing, this work establishes that the metamagnetic 
behavior of Kondo systems is intimately related to coher- 
ence. We have shown that the magnetization profile of 
a Kondo system is a useful probe of the onset of coher- 
ence at low temperatures. The location of the incipient 
plateau reveals the important energy scales; its height is 
a direct measure of the Luttinger sum rule. 
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APPENDIX A: BARE DENSITY OF STATES Hove singularity at w = 0: 



Numerical evaluation of the bare DOS is readily ac- 
complished using a power series expansion about the Van 



Po(w) 
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2ttH 



log 
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log 



9 / m 



25 / m 
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Truncating the expansion at eighth order muj/t results in 
a relative error of at most 10^^ in the region —t<uj<t 



and 10 ^ in the region — 4t < w < 4i. 
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